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ABSTRACT 

The present study the free convection in unsteady Couette flow of a viscous incompressible fluid confined between two vertical 
parallel plates in the presence of thermal radiation with heat source in the presence of uniform magnetic field is presented. The flow is 
induced by means of Couette motion and free convection currents occurring as a result of application of constant heat flux on the wall 
with a uniform vertical motion in its own plañe while constant temperature on the stationary wall. The fluid considered here is a gray, 
absorbing-emitting but non-scattering médium, and the Rosseland approximation is used to describe the radiative heat flux in the 
analysis. The dimensionless governing partial differential equations are solved by using regular perturbation technique. The results for 
the velocity, temperature and the skin-friction are shown graphically. The effects of different parameters are discussed. 
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1 INTRODUCTION 

Unsteady free convection flows in a porous médium have 
received much attention in recent time due to its wide 
applications in geothermal and oil reservoir engineering as 
well as other geophysical and astrophysical studies. Moreover, 
considerable interest has been shown in radiation interaction 
with convection for heat and mass transfer in fluids. This is 
due to the significant role of thermal radiation in the surface 
heat transfer when convection heat transfer is small, 
particularly in free convection problems involving absorbing- 
emitting fluids. The unsteady fluid flow past a moving píate in 
the presence of free convection and radiation were studied by 
Grief et al. [8] and Makinde [13]. From the previous literature 
survey about unsteady fluid flow, we observe that little papers 
were done in porous médium. The effect of radiation on MHD 
flow and heat transfer must be considered when high 
temperatures are reached. El-Hakiem [7] studied the unsteady 
MHD oscillatory flow on free convection-radiation through a 
porous médium with a vertical infinite surface that absorbs the 
fluid with a constant velocity. Cookey et al. [6] researched the 
influence of viscous dissipation and radiation on unsteady 
MHD free -convection flow past on infinite heated vertical 
píate in a porous médium with time-dependent suction. 

The fluid flow between parallel plates by means of Couette 
motion is a classical fluid mechanics problem that has 
applications in magnetohydrodynamic (MHD) power 

IJER@2013 



generators and pumps, accelerators, aerodynamics heating, 
electrostatic precipitation, polymer technology, petroleum 
industry, purification of crude oil, and also in many material 
processing applications such as extrusión, metal forming, 
continuous casting, wire and glass fiber drawing, etc. This 
problem has received considerable attention in the case of 
horizontal parallel plates, Attia, Choi et.al Kumar et. al and 
Soundalgekar [2, 4, 12, 16] than vertical parallel plates. 

In many engineering applications such as cooling of electronic 
equipments, design of passive solar systems for energy 
conversión, cooling of nuclear reactors, design of heat 
exchangers, chemical devices and process equipment, 
geothermal systems, and others. However, very few papers 
deal with free convection in Couette motion between vertical 
parallel plates. Jha [11] Fully-developed laminar free 
convection Couette flow between two vertical parallel plates 
with transverse sinusoidal injection of the fluid at the 
stationary píate and its corresponding removal by constant 
suction through the píate in uniform motion has been analyzed 
by Jain and Gupta [10]. In their study, the moving wall is 
thermally insulated and the wall at rest is kept at a uniform 
temperature. Narahari [14] Effects of thermal radiation and 
free convection currents on the unsteady couette flow between 
two vertical parallel plates with constant heat flux at one 
boundary. 
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The study of heat generation/absorption effects in moving 
fluids is important in view of several physical problems, such 
as fluids undergoing exothermic or endothermic chemical 
reactions. Hosssain et al. [9] studied the problem of natural 
convection flow along a vertical wavy surface with uniform 
surface temperature in the presence of heat generation / 
absorption. Alam et al. [1] studied the problem of free 
convection heat and mass transfer flow past an inclined semi- 
infinite heated surface of an electrically conducting and steady 
viscous incompressible fluid in the presence of a magnetic 
field and heat generation. Chamkha [3] investigated unsteady 
convective heat and mass transfer past a semi-infinite porous 
moving píate with heat absorption. 

The aim of the paper is to provide an exact analysis of 
unsteady free convection in Couette motion between two 
vertical parallel plates in the presence of thermal radiation in 
the presence of uniform magnetic field where the moving píate 
is subject to constant heat flux and the píate at rest is 
isothermal. These solutions are useful to gain a deeper 
knowledge of the underlying physical processes and it 
provides the possibility to get a benchmark for numerical 
solvers with reference to basic flow configurations. 
2 MATHEMATICAL ANALYSIS 

Consider the unsteady free -convective Couette flow of an 
incompressible viscous radiating fluid between two infinite 
vertical parallel plates in the presence of uniform magnetic 
field separated by a distance h . The x — axis is taken along 
one of the plates in the vertically upward direction and the 
y' — axis is taken normal to the píate. Initially, at time t' < , 
the two plates and the fluid are assumed to be at the same 
temperature T h and stationary. At timeí' > 0, the píate at y' 

= starts moving in its own plañe with an impulsive velocity 
U and is heated by supplying heat at constant rate whereas the 
píate at y' = h is stationary and maintained at a constant 

temperature T h . It is also assumed that the radiative heat flux 

in the x — direction is negligible as compared to that in the 
y — direction. As the plates are infinite in length, the 
velocity and temperature fields are functions of y and 

f only. Then under the usual Boussinesq's approximation, the 
flow of a radiating fluid is shown to be governed by the 
following system of equations: 

du' 



(„, ,\ du' crB () ¿ , v , 

— = g\T -71 \ + v — r — u u 

dt' \ ) dy' 2 p K 

(1) 
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dT ^T_5^ / A 
pCp dt'~ 8y' 2 dy' Qo V M 

(2) 

With the following boundary conditions 

t'<0: u' = 0,T' = T¡ t for 0<y'<h 



, , dT' q 

t >0: u =U, = — 

dy' k 



at y' = 
u' = 0, T = T¿ at y' = h 



(3) 

where g is the acceleration due to gravity, P the volumetric co- 
efficient of thermal expansión, V the kinematic viscosity, p 

the density, k the thermal conductivity, C the specific heat at 
constant pressure, q the constant heat flux, q r the radiative 
heat flux in y' — direction, T the fluid temperature, and u 
is the fluid velocity. 

The radiative heat flux term is simplified by making use of the 
Rosseland approximation Siegel and Howell [15] as 

4cr dT' 4 



9r 



(4) 



3k* dy' 



Where <7 is the Stefan - Boltzmann constant and k* is the 
mean absorption coefficient. It should be noted that by using 
the Rosseland approximation we limit our analysis to optically 
thick fluids. If temperature differences within the flow are 

sufficiently small such that T' 4 may be expressed as a linear 
function of the temperature, Then the Taylor series for T' 4 
about T h , after neglecting higher order terms, is given by 

T' 4 = 4T*T'-3T* 

(5) 

It is emphasized here that equation (5) is widely used in 
computational fluid dynamics involving radiation absorption 

problems Chung [5] in exprés sing the term T' 4 as a linear 
function. 

In view of Equations (4) and (5), Equation (2) reduces to 



(6) 



dT , d 2 T 16a dT! 3 d 2 T 



= k 



+ ■ 



dt' dy' 2 3k* dy' dy 



i2 



In order to solve the governing equations in dimensionless 
form, we introduce the following non-dimensional quantities: 
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U (hq/k) Ukv 



%-M = o 



(14); 0^-M=O 



(15) 



M = 



aB 2 h 2 K _h 2 ^_ Q h 2 
pv ' k pC p v 



,Pr = ^^,7? 



kk* 
4o-Tf 



(7) 

where Gr is the thermal Grashof number, Pr the Prandtl 

number, R the radiation parameter, <f) heat source parameter, t 

the dimensionless time, u the dimensionless velocity, y the 
dimensionless coordínate axis normal to the píate, \x the 
coefficient of viscosity and 8 is the dimensionless temperature. 
Then in view of Equations (7), Equations (1), (6) and (3) 
reduces to the following non-dimensional form of equations: 



d 2 u du 



1 



dy 2 dt 
(8) 



Mu u + Grd 

K 



d 2 e 



de 



(37? + 4) — -37? Pr 37? Vrád 

V ' dy 2 dt y 

(9) 

The initial and boundary conditions are 
t<0: u = O,0 = O for 0<y<l 

t>0: u = l, — = -1 

dy 



m = 0,0 = at 



at y - 
y = l 



(10) 



3 SOLUTION OF THE PROBLEM 

Equation (8) - (9) are coupled, non - linear partial differential 
equations and these cannot be solved in closed - form using 
the initial and boundary conditions (10). However, these 
equations can be reduced to a set of ordinary differential 
equations, which can be solved analytically. This can be done 
by representing the velocity, temperature and concentration of 
the fluid in the neighbourhood of the fluid in the 
neighbourhood of the píate as 



f = f (y) + se n 'f l (y) +0(e 2 ) 



(11) 



Substituting (11) in Equation (8) - (9) and equating the 
harmonic and non - harmonic terms, and neglecting the higher 

order terms of O^í -2 ) , we obtain 

u" ) -j3 l u Q =-Gr0 o (12); ^"-/V^-M 
(13) 



The corresponding boundary conditions can be written as 



u n =l,u, — 0, — - = -l, — 1 - = 0, aty = 
1 8y 8y (16) 

u Q ->0,U! ->0, O ->0,^ ->0, asy^l 

The solutions of Equations (12) - (15) under the initial and 
boundary conditions (16) by perturbation technique is given 
by 

u(y,t) =Z,e" hy +Z 2 e" v +Z 3 e n%y +Z 4 e" hy 

0(y,t) = D l e" hy +D 2 e miy 
APPENDIX 

The constant not given Brevity of the space 
4 RESULTS AND DISCUSSION 

An exact solution to the problem of natural convection in 
unsteady Couette flow between two long vertical parallel 
plates in the presence of constant heat flux and thermal 
radiation have been presented in the preceding section. In 
order to get the physical insight into the problem, the 
numerical valúes of the temperature field, the velocity field, 
the skin-friction, the Nusselt number, the volume flow rate and 
the vertical heat flux are computed for different valúes of the 
system parameters such as Radiation parameter (R), Grashof 

number (Gr) , Prandtl number (Pr) and heat source 
parameter (^) . Figure (1) presents the velocity profiles for 
both air and water (Pr = 7.0) in the case of puré convection 
(R -> 100) for different valúes of Gr . It is seen that the 

velocity of air and water increases with increasing Gr . At a 
smaller Gr the velocity distribution is monotonic, but at a 
higher time it passes through a máximum near the moving 
píate when the buoyancy effect partly suppresses the inertial 
effects of the píate velocity. Moreover, the velocity of air is 
greater than the velocity of water. Physically this is possible 
because fluids with high Prandtl number have greater 
viscosity, which makes the fluid thick and henee move slowly. 
Figure (2) presents the velocity profiles of air for different 
valúes ofK . It is observed that the velocity increases with 
increasing K . Physically this is possible because as the 
Grashof number or time increases, the contribution from the 
buoyancy forcé near the moving hot píate become more 
significant and henee a small rise in the fluid velocity near the 
píate is observed. Figure (3) illustrate the influences of M 
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(Magnetic parameter) on velocity profiles respectively. It is 
found that the velocity decreases with increase of magnetic 
parameter M. It is also found that the velocity decreases away 
from the píate and becomes minimum and finally takes 
asymptotic valué. Finally here we also see that point of 
separation takes place for different valúes of magnetic 
parameter. Figure (4) present typical profiles for the velocity, 

for various valúes of a heat source respectively. As 
shown, the velocity decreasing with increasing (j) . In the 
event that the strength of the presence of a heat source 

effect causes a reduction in the thermal state of the fluid, thus 
producing lower thermal boundary layers. The effect of the 
Prandtl number on the velocity shown in figure (5). As the 
Prandtl number increases, the velocity decreases. Further, 
figure (6) it is observed that the fluid velocity decreases with 
increasing valué of R. This result may be explained by the fact 



that an increase in the radiation parameter R — 



kk* 

4< y 



for 



fixed k and T h means an increase in the Rosseland mean 

absorption coefficient k * . When radiation is present, the 
momentum boundary layer was found to be thicken, which is 
in agreement with the observation made earlier with regard to 
the temperature variations of air. Figure (7) present typical 
profiles for the temperature for various valúes of a heat source 

respectively. As shown, the temperature decreasing with 

increasing. Figure (8) shows that the temperature profiles for 

different valúes of Prandtl number (Pr) ■ It is observed that an 

increase in the Prandtl number results a decrease of the 
thermal boundary layer thickness and in general lower average 
temperature within the boundary layer. The reason is that 
increasing valúes of Prandtl number equivalent to increase the 
thermal conductivities and therefore heat is able to diffuse 
away from the heated píate more rapidly. Henee in the case of 
increasing Prandtl numbers, the boundary layer is thinner and 
the heat transfer is reduced. Boundary layer suction is the 
technique in which air pumps is used to extract the boundary 
layer at the wing. Further, increment of suction parameter 
decreases the fluid temperature. Figure (9) shows the 
temperature profiles for different valúes of the Radiation 
parameter R, it is noticed that an increase in the radiation 
parameter results decrease in the temperature with in boundary 
layer, as well as decreased the thickness of the temperature 
boundary layers. Figure (10) presents the skin-friction 



variation with M in the puré convection case for different 
valúes of Gr at the moving píate. It is observed that the skin- 
friction increases with increasing Gr . 
5 CONCLUSIONS 

The temperature of the fluid increases with increasing time 
whereas it decreases due to an increase in the valué of 
radiation parameter. In the case of puré convection (i.e. in the 
absence of radiation), the velocity of the fluid increases with 
increasing Grashof number, but falls owing to an increase in 
the Prandtl number. The velocity of the fluid increases with 
increasing Grashof number and time but it decreases owing to 
an increase in the valué of the radiation parameter. The skin- 
friction at the moving píate increases with increasing valúes of 
Grashof number and time for air flows. 

We may conclude therefore, that the interaction between 
the radiation, buoyaney forces and the applied shear induced 
by a uniform vertical motion of the hot wall can affect the 
configuration of the flow field significantly. 
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Figure 1 : Velocity profiles for different valúes of Gr (Puré convection case) 
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Figure 3: Velocity profiles for different valúes of M 
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Figure 4: Velocity profiles for different valúes of <| 
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Figure 5: Velocity profiles for different valúes of Pr 
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Figure 6: Velocity profiles for different valúes of R 
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Figure 2: Velocity profiles for different valúes of K 
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